Substructure Boosts to Dark Matter Annihilation from Sommerfeld Enhancement 
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The recently introduced Sommerfeld enhancement of the dark matter annihilation cross section 
has important implications for the detection of dark matter annihilation in subhalos in the Galactic 
halo. In addition to the boost to the dark matter annihilation cross section from the high densities 
of these subhalos with respect to the main halo, an additional boost caused by the Sommerfeld 
enhancement results from the fact that they are kinematically colder than the Galactic halo. If 
we further believe the generic prediction of CDM that in each subhalo there is an abundance of 
substructure which is approximately self-similar to that of the Galactic halo, then I show that 
additional boosts coming from the density enhancements of these small substructures and their 
small velocity dispersions enhance the dark matter annihilation cross section even further. I find 
that very large boost factors (10^ to 10^) are obtained in a large class of models. The implications 
of these boost factors for the detection of dark matter annihilation from dwarf Spheroidal galaxies 
in the Galactic halo are such that, generically, they outshine the background gamma-ray flux and 
are detectable by the Fermi Gamma-ray Space Telescope. 



I. INTRODUCTION 

In the standard flat, Gaussian, adiabatic, and scale- 
invariant ACDM paradigm ~ 80 % of the matter density 
of the Universe is in the form of dark matter (DM) [1]. 
The properties of dark matter are largely unconstrained, 
although a weakly interacting massive particle (WIMP) 
with a mass ~ 100 GeV - 1 TeV has many attractive 
features, including that it leads to a relic density of dark 
matter that is remarkably close to the measured value 
[2]. No interaction of dark matter either with itself or 
with standard model (SM) particles other than gravity 
has ever been observed. However, at the present day we 
cannot exclude that the dark sector is brimming with a 
rich phenomenology at weak-interaction scales or lower. 

One exciting possibility is that dark matter self- 
annihilates at a level observable in the Galactic halo 
today. Detecting this self-annihilation seems unlikely 
considering the standard self-annihilation cross section; 
(c'"'^)ann ^ 3 X 10~^^ cm^ s~^ for consistency with the 
observed relic density. However, there are various mech- 
anisms to boost the annihilation signal to a value larger 
by a few orders of magnitude, e.g., dumpiness of the 
dark halo [3, 4, 5, 6] or a combination of coannihila- 
tions and Sommerfeld enhancement [7, 8]. The possi- 
bility of detecting the DM annihilation signal has at- 
tracted much attention recently because of the anoma- 
lous excesses of high energy electrons and positrons re- 
ported by two different experimental groups: the Pay- 
load for Antimatter Exploration and Light-nuclei Astro- 
physics (PAMELA) satellite, which saw a sharp upturn 
in the positron fraction e+/(e+ + e~) from 10 to 100 
GeV [9], and the Advanced Thin Ionization Calorime- 
ter (ATIC) group, which reported an excess over the 
expected background in the combined number of elec- 
trons and positrons at energies 500 to 800 GeV [10]. The 
PAMELA result on the positron fraction confirms earlier 
results [11, 12, 13, 14, 15, 16, 17] and can be explained 
by DM annihilation [18, 19, 20]. There is still a good 



chance that these anomalies can be obviated by conven- 
tional astrophysical effects, e.g., pulsars [21, 22, 23, 24]. 
Nevertheless, a plethora of models has been proposed in 
the last few months explaining the PAMELA and ATIC 
results in terms of DM annihilation into electron-positron 
pairs (e.g., [25, 26, 27, 28, 29, 30, 31]), a description that 
also naturally accounts for other experimental anomalies, 
such as the "WMAP haze" [32, 33, 34], EGRET obser- 
vations of an excess of high-energy gamma-rays in the 
Galactic center [35], and recent HESS measurements of 
the flux of very high-energy electrons [36]. Model build- 
ing is complicated by the fact that PAMELA did not 
observe a similar excess in the ratio of anti-protons to 
protons [37], which means that the coupling of the DM 
particles to quarks must somehow be suppressed. 

The largest obstacle to modeling the PAMELA/ ATIC 
signal with DM annihilation is the very large boost fac- 
tor (10^ to 10^) with respect to the fiducial (ct'u) ann ^ 
3 X 10~^^ cm^ s~^ required to produce an observable 
annihilation flux. The boost factor from the dumpi- 
ness of the dark matter halo is expected to be small 
[5, 6, 38, 39], typically of the order 1 to 10. Therefore, 
a different mechanism must be responsible for the large 
boosts necessary. One such mechanism that has received 
much attention lately is the Sommerfeld enhancement 
[8, 40, 41, 42, 43, 44, 45, 46, 47, 48]. The Sommerfeld 
enhancement is a generic effect present whenever there is 
an attractive force acting between the dark matter par- 
ticles, e.g., a Yukawa interaction or a gauge interaction 
through vector bosons [26]. The main effect of this at- 
tractive force is to enhance the annihilation cross section 
(TV with a factor proportional to I3~'^ {P = v/c); this is 
known as "1/v" enhancement. There are regions in the 
parameter space made up of the dark matter mass m^, 
the force carrier mass m^, and the force "fine-structure 
constant" a = coupling^ /47r in which the enhancement 
is near a resonance, where it approximately behaves as 
However, the enhancement saturates at a certain p 
because of the finite range of the attractive force — large 
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{a ^ 10~^ to 10~^) attractive forces with infinite range 
lead to a burst of DM annihilation in the first DM halos 
formed at z ~ 100 - 200 which would lead to observable 
effects incompatible with measurements of the diffuse ex- 
tragalactic gamma-ray background today and the cosmic 
microwave background [49]. Since the typical velocities 
of dark matter particles in the Galactic halo today are 
of the order of the velocity dispersion of the halo, which 
is non-relativistic {P ~ 10~^), the cross section of dark 
matter particles annihilating in the halo could be signif- 
icantly enhanced by this Sommerfeld effect. 

Both observations [50, 51, 52, 53, 54] as well as nu- 
merical simulations [38, 55, 56, 57, 58, 59, 60, 61] have 
firmly established that dark matter halos in CDM are 
not smooth structures, but that they are filled with a 
large amount of substructure. These substructures are 
both in the form of subhalos, self-bound objects orbit- 
ing in the Galactic halo, and in the form of streams, 
relics of tidally disrupted subhalos. Since these sub- 
structures are over dense regions in the halo, the dark 
matter annihilation cross section will be correspondingly 
higher: Annihilation cross sections scale as the density 
squared. The effects of these kinds of substructure boosts 
have been studied extensively during the last decade 
[5, 6, 39, 62, 63, 64, 65, 66, 67, 68]. However, these subha- 
los are not only denser than the Galactic halo in which 
they orbit, they are also kinematically colder than the 
Galactic halo, with velocity dispersions that are typically 
two orders of magnitude smaller [69, 70]. This means 
that the Sommerfeld enhancement in these subhalos can 
be many orders of magnitude larger than the Sommer- 
feld enhancement acting in the Galactic halo taken as a 
whole. 

In addition to this larger boost caused by the small 
velocity dispersion in the subhalos, these subhalos them- 
selves are filled with an abundance of substructure which 
is approximately self-similar to that of the main halo. 
Each of these even smaller substructure has a larger den- 
sity than the subhalo in which they are orbiting (quanti- 
fied by the well established concentration-mass relation- 
ships for CDM halos [71, 72]) as well as a smaller velocity 
dispersion. This opens up the possibility that the dark 
matter annihilation signal in these already cold subha- 
los can be boosted even further by contribution of their 
substructure. 

In this paper we quantify the magnitude of this new 
substructure boost and discuss its implications for the 
detection of dark matter annihilation from subhalos in 
the Galactic halos. As we will show, the total combined 
boost from Sommerfeld enhancement and substructure 
can reach very large values (10^ to 10^), and these boosts 
are only weakly dependent on the assumptions that we 
have to make about the abundance of substructure and 
the DM density profile. We use these total boosts to- 
gether with published DM density profiles of some of the 
dwarf Spheroidal satellite galaxies of the Milky Way to 
predict the total flux of gamma-rays from DM annihi- 
lation that would be detectable by the Fermi Gamma- 



ray Space Telescope (Fermi) and find that in most of 
the models that posit the Sommerfeld enhancement the 
DM annihilation signal from the dwarf Spheroidals shines 
above the diffuse extragalactic gamma-ray background, 
which is the only background signal at the high Galac- 
tic latitudes of these satellite galaxies. Finally, we also 
consider the prospects for detecting smaller subhalos of 
the Galaxy by computing annihilation fluxes of subhalos 
in the Via Lactea II simulation. Whether these subhalos 
will be detectable in the near future depends on the exact 
details of the Sommerfeld enhancement. 



II. THE SOMMERFELD ENHANCEMENT 

The simplest example of the Sommerfeld effect that 
contains all of the generic features described in the in- 
troduction is that from an attractive Yukawa type force 
mediated by a spin-0 boson (j) between the DM particles, 
and we adopt that description here. The Sommerfeld ef- 
fect is a non-perturbative effect in the sense that in the 
language of quantum field theory it comes about from 
Feynman diagrams in which the force carrier is exchanged 
many times before the annihilation actually happens, giv- 
ing rise to so called "ladder diagrams" . However, we can 
study the Sommerfeld effect in simple, non-relativistic 
quantum mechanics (see the appendix in Ref. [26]). We 
are instructed to solve the radial Schrodinger equation 
(in natural units) 



1 d2^(r) 



rriy 



dr^ 



+ _e-^^^^(r) = -m^(3^ilj{r) , (1) 



subject to the boundary condition di/j/dr = im^f3ilj as 
r ^ oo. Here and are the mass of the dark matter 
and the force carrier, respectively. Using the substitution 



r am^r this becomes 



d^?/^(r) 



dr^ 

with boundary condition 
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The Sommerfeld enhancement S is then given by 



(2) 



(3) 



(4) 



Equation (2) shows that the Sommerfeld enhancement 
is only a function of two variables, which can be chosen 
as a/^ and am^/m^. The qualitative features of equa- 
tion (2) have been widely discussed and, therefore, we 
will merely mention them here without proof (see, e.g., 
Refs. [8, 26]). When the kinetic energy of the collision is 
much larger than the product of the force carrier mass 
and the coupling, then the Yukawa interaction can be 
approximated as an attractive Coulomb interaction, for 
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which an analytical solution in terms of hypergeometric 
functions exists. The Sommerfeld enhancement is simply 
given by 



-T^a/ (3 



(5) 



From this we see why the Sommerfeld enhancement is 
often referred to as a "1/'^" enhancement for small ve- 
locities /3 <C 1. In the limit of a vanishing interaction 
{a/ j3 <C 1) the Sommerfeld enhancement disappears as 
5* ^ 1. The 1/v behavior, which diverges for small (3 
and leads to an overproduction of high-energy photons 
in the early Universe [49], saturates as a consequence of 
the finite range of the Yukawa interaction at a velocity 
P ^ m(f)/m^ and settles on a value S ^ am^/m^f). How- 
ever, at some parameter combinations the Yukawa po- 
tential has bound states, at which there are large boosts 
to the annihilation cross section which go approximately 
as 5 oc l3~'^. These resonance regions, however, are cut- 
off at small P by the finite lifetime of the bound state. 
We expect these three different behaviors {l/v, satura- 
tion at low V, and resonance effects from the existence 
of bound states) to be generic among different models of 
the Sommerfeld enhancement (see also the discussion in 
the appendix of Ref . [26] ) . 

Equation (2) can also be numerically integrated effi- 
ciently, by integrating the solution ?/^(r) = e*^^/" inwards 
from a radius satisfying 



1 



v2 • 



(6) 



In the Coulomb regime {m^f]'^ ^ OLTricj)) this means that 
the initial value for the integration should be such that 
r ^ o? I otherwise we should start at r ^ am^/rricf). 
Figure 1 shows the Sommerfeld enhancement in the pa- 
rameter plane made of a and mcjy/m^. All of the qual- 
itative features we discussed above are present: In the 
Coulomb regime (lower part of the two panels) we see 
that the enhancement scales as a, saturating at 5* ~ 1 
when a becomes small; The fact that the upper left cor- 
ner of both panels is approximately the same shows that 
the Sommerfeld enhancement saturates at ~ m^/m^; 
The sharp diagonal lines in the upper right of both pan- 
els show the resonance region, in which the enhancement 
is sharply peaked. Because of the scaling properties of 
equation (2), shifting log^^ is equivalent to shifting both 
log a and log m^/m^ by the same amount. Therefore, 
the (5 behavior of the enhancement can also be read off 
from these figures along diagonals. Inspection of the res- 
onance diagonals in the second panel of Figure 1 shows 
that the enhancement grows as in these regions. 

Of course, when considering the Sommerfeld enhance- 
ment of a halo of a certain velocity dispersion cjy we need 
to average the Sommerfeld enhancement over the distri- 
bution of relative velocities of that halo. We will approx- 
imate the one-particle velocity distribution of the halo as 
a single truncated Maxwell-Boltzmann distribution with 
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FIG. 1: Sommerfeld enhancement as a function of the parti- 
cle physics parameters for two different values of the relative 
velocity [3. Because of the scaling properties of equation (2) 
the behavior of the Sommerfeld enhancement as a function of 
log^Q [3 can actually be read off along diagonals in the plane. 
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(7) 



This approximation of a single Maxwell-Boltzmann dis- 
tribution for the velocity distribution means that we as- 
sume that the halo has a constant velocity dispersion over 
the halo, i.e., that the halo is isothermal. This is not nec- 
essarily the case — indeed, it is not if we believe that the 
density profiles of DM halos are well fit by Einasto or 
NFW profiles — but it is not a bad assumption for inter- 
mediate distances from the halo's center, i.e., between 
the innermost and outermost parts of the halo [73]. A 
more detailed treatment of the Sommerfeld enhancement 
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in subhalos should take the non-uniform velocity disper- 
sion in subhalos into account. Note that in the pure \ jv 
regime the averaging will tend to sustain the Xjv behav- 
ior up to a factor as we have that 



(8) 

since the truncation does not significantly infiuence the 
result of these integrals. Similarly, in the resonance re- 
gion the \jv^ behavior is conserved by the averaging. 
When the enhancement has leveled off the averaging has 
no effect. 

In Figure 2 we show this averaged Sommerfeld en- 
hancement for four different decades of velocity disper- 
sion. The overall behavior of the Sommerfeld enhance- 
ment seen in Figure 1 is indeed conserved well by the 
averaging: this has to do with the fact that the most 
probable relative velocity in a Maxwell-Boltzmann rela- 
tive velocity distribution is the velocity dispersion. The 
Sommerfeld enhancement attains very large values for 
subhalos with very small velocity dispersions, mostly, but 
not exclusively, in the resonance regions. This means 
that the smallest subhalos that exist in CDM will get 
large boosts to their DM annihilation from Sommerfeld 
enhancement. 



III. SUBSTRUCTURE BOOSTS FROM 
SOMMERFELD ENHANCEMENT 

The photon flux from DM annihilation from a solid 
angle lS£l along a given line-of-sight (los) is given by 



dTV^ 
A A At 



dN^^i {av)i 



dE 



S {m^/m^,a,ay) 



X / -r I (0: 

I An 47r Jios 



(9) 



where £^th is the threshold energy of the detector, 
i denotes (possibly) different final states of the DM 
annihilation, is the spectrum of the annihila- 

tion to that state, {av)i is the annihilation cross sec- 
tion to state i (without the Sommerfeld enhancement), 
S {mcjy/m^^a, (3) is the Sommerfeld enhancement which 
we have factored out of the cross section, and p{l) is the 
DM density along the line-of-sight. We emphasize that 
this framework also allows us to consider different final 
annihilation products, such as the electrons and positrons 
observed by PAMELA/ATIC. 

We will focus our analysis on the "structure quantity" 
£(M), which we define here as 

C{M)=S{m^/m^,a,a,) f dl p\l) (10) 

This structure quantity does not only depend on the mass 
and internal properties of the DM halos in which we are 



interested, it depends on the particle physics parameters 
as well. We will keep all of the other factors in equation 
(9) fixed. However, in order to derive actual numbers for 
the photon fiux from DM annihilation, we briefiy describe 
the other quantities appearing in equation (9). 



A. Some details of the particle physics model 

For the annihilation cross section prior to any Som- 
merfeld enhancement we will adopt the fiducial value 
{cFv) = 3 X 10~^^ cm^ s~^, and we choose a fiducial mass 
for the DM particle of = 700 GeV. The energy spec- 
trum of the DM annihilation depends on the model used 
to describe the DM and its interactions; it is outside of 
the scope of this paper to provide a detailed treatment 
of this for different DM models. Therefore, we will adopt 
an annihilation spectrum for DM annihilation into pho- 
tons appropriate for the newer models that were proposed 
to explain the PAMELA/ATIC results. Specifically, we 
use a gamma ray spectrum from a cascade annihilation 
model. In these models the DM annihilates into a new 
light degree of freedom which then in turn decays into 
leptons and photons. The injection spectrum per DM 
annihilation is given by [74] 



dN^ _a^M 1 + (1 - E/m^y 



dE 



E 
In (4(1 



^/m,))-21n(^ 



(11) 



in which is the electromagnetic coupling constant 

and mi is the mass of the leptons created in the process 
XX ~^ Integrating this spectrum using an energy 

range appropriate for Fermi we find that 



250 GeV 



4 GeV 



dE 



0.5. 



(12) 



We find a similar number if we integrate the spectrum 
from one cascade step. It is interesting to note that these 
photon spectra have fewer photons in them, about an or- 
der of magnitude less, than the annihilation spectra from 
neutralino annihilation which are often used [75, 76]. Of 
course, since the energy spectrum contributes a multi- 
plicative factor to the DM annihilation fiux, the effect 
of different annihilation spectra can be incorporated by 
simply multiplying our results with the appropriate fac- 
tor. 



B. The Dark Matter profile 

Since the photon fiux from DM annihilations depends 
on the density squared of the DM the exact form of the 
inner density profile has important ramifications for the 
detectability of DM annihilation. The density profile of 
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dark matter halos has been discussed extensively in the 
last few years. While the shape of the outer density 
profile seems to have converged in the various numeri- 
cal simulations of DM halos and a consensus has been 
reached, the detailed shape of the inner density profile is 
still a matter of debate. The main point of contention is 
whether the density profile is cusped near the center of 
the DM halo, and if so, how steep this cusp is, or whether 
the density asymptotes to a constant value at the center. 
Since this discussion has a non-negligible impact on the 
prospects for DM annihilation detection, we will briefly 
discuss the various plausible density profiles. 

The inner structure of DM halos is one of the key pre- 
dictions of the CDM paradigm, and over a decade ago it 
was shown that the spherically averaged density profiles 
of DM halos are well described by a universal profile, the 
so-called NFW profile [77, 78, 79], depending only on the 
virial mass of the halo and its concentration. Moreover, 
the concentration was shown to depend systematically 
on the halo mass [71, 72, 78, 79, 80, 81]. This profile 



was found to provide a good fit to many numerical sim- 
ulations (e.g., [82, 83, 84]). While there seems to be 
a consensus between different numerical simulations on 
the shape of the outer density profile, no agreement was 
reached on the inner density structure of the DM halos 
found in the simulations. The NFW profile has an inner 
slope of -1, other studies, however, found evidence for 
steeper inner slopes [85, 86, 87, 88, 89, 90, 91], shallower 
inner slopes [92], or didn't find any convergence toward 
a power-law behavior [93]. Even the universality of the 
density profile has been brought into question [94, 95]. 

More recently, high-resolution numerical simulations 
have essentially ruled out the steepest cusps, such as 
Moore's profile, which has an inner cusp a with logarith- 
mic slope of -1.5, however steeper inner cusps than the 
original NFW profile are still observed in numerical sim- 
ulations [96]. Other recent simulations find inner slopes 
that are shallower than an NFW profile [97, 98, 99]. It 
also has become apparent in recent years that the inner 
density profile of DM halos might not be described by a 
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power-law divergence, but that, on the contrary, the log- 
arithmic slope becomes ever shallower, leading to a cored 
density profile [61, 98, 99, 100, 101, 102, 103, 104]. There- 
fore, density profiles such as the Einasto profile [105, 106] 
have been shown to provide an excellent fit to the numer- 
ical simulations with the highest resolution yet [61, 98]. 

The implications of the exact form of the inner density 
profile have been widely discussed before in the context 
of various DM annihilation models (e.g., [75, 107, 108, 
109, 110, 111, 112, 113, 114, 115, 116]). Roughly speak- 
ing, the annihilation flux is smaller for inner density pro- 
files which are cored, and larger for inner density profiles 
which show a cuspy behavior. The steeper the slope of 
the inner cusp, the larger the DM annihilation flux (for 
obvious reasons). Since the core vs. cusp debate is still 
undecided we will consider DM distribution models of 
both types. 

The cuspy DM profile we will use is a generalized NFW 
profile (GNFW): 



p{r) 



(1 + rf~ 



r/vs 



(13) 



in which ps is a characteristic density and a charac- 
teristic (scale-)radius. The original NFW profile is re- 
covered for 7 = 1 [79], which we will use as the fiducial 
model in what follows. Varying the value of the inner 
logarithmic slope 7 will allow us to describe the effect of 
steeper and shallower inner slopes than the NFW model. 
Note that this functional form does not include the exact 
form of Moore's profile; this need not concern us since in- 
ner cusps as steep as Moore's profile are excluded by the 
latest numerical simulations. 

The contribution of the DM density profile on the 
structure quantity C can be simply evaluated in the con- 
text of any density model. As in Ref. [65] we assume 
that the distance to the subhalo D is much larger than 
the scale-radius of the DM halo. For an NFW profile 
90% of the flux originates from the region within one 
scale-radius, and we will, conservatively, only take con- 
tributions to the DM annihilation flux coming from this 
region for all of the different profiles. Therefore, we find 
that 



An 47r 



los 



(14) 



The normalization of this relationship depends on the 
logarithmic slope of the inner cusp and is given by 



Afi{x) = 47r ■ 



2x-5 (_16+llx-2x2) 

-30 + 47 X - 24^2 + 4^3 



(15) 



The normalization for the NFW profile is 77r/6 and it 
depends strongly on the value of the logarithmic slope 
of the inner density profile, as can be seen in Figure 3. 
The bottom panel of Fig. 3 shows the fraction of the DM 
annihilation flux coming from region within one scale- 
radius of the halo. For reasonable values of the inner 
slope this fraction is in the range 80 to 95%. 




FIG. 3: Top: Normalization of the C oc p^r^/D^ relation of 
equation (14) for the GNFW profile. The functional depen- 
dence of the normalization on the exponent 7 in the profile 
given in equation (13) is given in equation (15). Bottom: 
Fraction of the emission that comes from the region within 



The Einasto profile is characterized by an ever- 
decreasing logarithmic inner slope. It's functional de- 
pendence on the radius is given by 



p{r) = p-2 exp 



(16) 



Here, p-2 is the density at the radius r_2 where the local 
slope is -2 and a is a shape parameter often fixed at a 
value ~0.17 [61]. As above, assuming D r_2, calcu- 
lating the flux originating from within the radius r_2, we 
find that 
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Z)2 



(17) 



The normalization depends on the value of a and is given 
by 



7 



3 4 

a a 



(18) 



where 7[a,^] is the lower incomplete gamma function. 
The dependence of this normalization on a is shown in 
Figure 4. In order to compare the Einasto profile as de- 
scribed here to the GNFW profile, we can relate the pa- 
rameters (ps, Ts) of the GNFW profile to the parameters 
(p_2,r_2) of the Einasto model by calculating the ra- 
dius at which the logarithmic slope of the density versus 
radius profile equals -2. Thus, we find that 



r-2 Ts 



P-2 



23-7 



(19) 
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FIG. 4: Top: Normalization of the C oc p^_2^^-2/D^ relation 
of equation (17) for the Einasto profile. The functional de- 
pendence of the normalization on a in the profile given in 
equation (16) is given in equation (18). Bottom: Fraction of 
the emission that comes from the region within r-2- 



Note that the DM annihilation flux from the Einasto 
profile, although it seems larger from Figs. 3 and 4, is 
actually smaller in general, since there is a factor of 2^~^^ 
involved in going from one to the other. 

Further on it will sometimes be useful to describe the 
DM density profiles by different parameters than a char- 
acteristic density and radius. Another set of parameters 
is given by a concentration parameter and a mass. The 
virial radius ryir of a halo is defined to be the radius 
within which the mean density of the halo is A times the 
critical density pc, where A = 1781^^'^^ ^ 100 [117, 118]. 
The virial mass is simply the mass contained within a 
sphere of radius Tyir 
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(20) 



The concentration for both the GNFW profile as for the 
Einasto profile is then defined by 



r-2 



(21) 



A characteristic mass of a halo is given by the virial mass 
defined above. The virial mass of the GNFW profile in 
terms of the characteristic density ps and the character- 
istic radius r. is 



(22) 



in which the function /(c) is given by 

.3-7 



3-7 



2Fi(3-7,3-7;4-7;-c) . (23) 



The function 2F1 (a, 6; c; z) is the Gauss hypergeometric 
function. In the 7 = 1 case /(c) is given by the familiar 
form /(c) = ln(l + c) — c/(l + c). The virial mass for the 
Einasto profile in terms of the characteristic density p-2 
and the characteristic radius r_2 is 
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31na + 2-ln^ 



a 



7 



-,-c^ 
a a 



(24) 

In this last expression 7 [a, is the lower incomplete 
gamma function. The density can then be expressed in 
terms of the concentration by using the definition of the 
virial mass given in equation (20). 

Another pair of parameters that one can use to char- 
acterize the DM density profile is the maximum circular 
velocity 'i^max and the radius at which this maximum oc- 
curs, Tmax- These parameters are especially well-suited 
when dealing with numerical simulations, since they can 
be readily read off from the raw simulation data. The 
circular velocity is given by 



(25) 



in which x = r/r_2 and f{x) is defined in equation (23) 
in the case of a GNFW profile and is given by 7 ^x"] 
in the case of an Einasto profile. One can find rmax by 
(numerically) solving the equation 



_d_ 

dx 



0, 



(26) 



and 'Umax follows then immediately. 



C. Calculating the total boost 

As is customary, we will write the total structure quan- 
tity jC{M) as a sum of a smooth component and a sub- 
structure component: 



C{M) = Csm{M) + C,^^{M) . 



(27) 



The smooth contribution to the structure quantity is the 
product of the Sommerfeld enhancement S and the line- 
of-sight integral over the DM density squared, given for 
the different DM density profiles in equations (14) and 
(17). Therefore, 



{M)=S 



(28) 



in which Csm,v=c is the luminosity that would come from 
the smooth component of the halo in the absence of Som- 
merfeld enhancement. 
We can now write 



£(M) = SCsra,v=c{M) + 4ub(A^) 
Aub(M) 
Csm,v=c{M) 



(29) 
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We will from now on write the substructure boost as 
B{M). Then we have the following 



B{M) 



1 / . dTV./ N 

-— - / dm - — £{m) 

y^c{M) J dm ^ ^ 



(M) Csm,v^c{M) 

-—^ / dm — [^(m) + 5(m)] £sm,^=c(^) . 

f^^m.,v—c J CLm 

(30) 



Here we have introduced the subhalo mass function 
dA^/dm, which gives the abundance of subhalos of mass 
m. 

Before proceeding we should specify the upper and 
lower limits of the integral appearing in equation (30). 
The integration over ever smaller substructures is cut- 
off at the size of the smallest subhalos. The mass 
scale at which this cut-off occurs is the thermal free- 
streaming scale, a DM model-dependent quantity, which 
in many DM models lies in the range 10~^ to 10~^^ 
[119, 120, 121, 122, 123]. The upper limit to this integral 
is some fraction q of the halo mass M, since the abun- 
dance of substructure does not continue up to this mass 
scale. In the following we will set q equal to 0.1. 

In order to calculate the total boost from substructure 
to the annihilation cross section, we will write all of the 
quantities appearing in equation (30) as a function of 
the mass m. For instance, from numerical simulations, 
the subhalo mass function is well constrained to follow a 
power-law type behavior: 



dA^ 
dlnm 



(X 



/ m \ ^ 



(31) 



The value of the exponent n is still under debate. 
The two highest resolution A'-body simulations yet, the 
Aquarius and Via Lactea projects, both find values for 
n that put roughly equal amounts of mass per subhalo 
mass decade. While the Via Lactea simulation seems to 
prefer a value n = — 1 [6, 124], which would lead to a log- 
arithmically diverging total mass in substructures if the 
subhalo abundance were not cut off at the thermal free- 
streaming limit of the dark matter, the Aquarius project 
prefers a value n = —0.9 [61]. As the Aquarius project 
has made detailed convergence studies that support their 
claim for a value of n that is shallower than -1, we will 
adopt the value n = —0.9 as our fiducial value, although 
even shallower slopes have also been observed [125]. The 
same slope has also been observed in other simulations 
[126]. 

With the exponent n fixed, we turn to the normal- 
ization of the relation given in equation (31). Again, 
different high-resolution numerical simulations give dif- 
ferent answers to this question as it relates to the abun- 
dance of subhalos of Galactic subhalos. The problem 
of the normalization is complicated by the fact that 
numerical simulations have only recently begun to be 
able to resolve substructure in subhalos in a Galaxy-size 
halo. While the simulations generally agree that the re- 
solved substructure abundance of subhalos in the host 



halo adds up to ~ 10% of the total mass of the halo 
[55, 58, 64, 91, 127, 128], the question of whether the 
substructure abundance in subhalos is just a scaled-down 
version of the substructure abundance in the host halo re- 
mains unresolved. One would not expect the abundance 
to be self-similar, since, while substructure in both the 
main host halo and subhalos gets diminished by tidal 
disruptions, substructure in the subhalos does not get 
replenished by the infall of halos from the field, as is the 
case for the host halo. Nevertheless, many numerical sim- 
ulations do find approximately self-similar substructure 
abundances [57, 129]. Recently, the Aquarius simula- 
tion has concluded from its simulation that the substruc- 
ture abundance is not self-similar, although the deviation 
from the self-similar relation is small [61]. Therefore, we 
will adopt a normalization here which places 10% of the 
mass of the subhalo in the form of "resolved" substruc- 
ture, where we adopt the definition that resolved sub- 
structure is in halos of mass > 10~^M. 

In the previous section we calculated £sm,v=c foi" the 
different DM density profiles and found in all cases that 



,(m) = 



,(p(m),r(m)) , 



(32) 



for a characteristic density p and radius r, the exact defi- 
nition of which depends on the specific DM density profile 
used. In order to find the mass dependence of >Csm,v=c7 
we use the description of the DM density profile in terms 
of virial mass M^\^ and concentration c. We then make 
use of the fact that many simulations have found a rela- 
tion between the concentration of halos and their mass 
[71, 72, 80, 81], such that everything depends only on 
the mass. This relation is well established over a large 
range in halo masses, and we will adopt the particular 
relationship found in Ref. [130] 



10.5 



1012 Mo 



-0.11 



(33) 



in good agreement with the results of other, similar 
analyses [71, 131], although the normalization of the 
concentration-mass relation is about 15% lower than the 
original relation. However, the concentration of subha- 
los is generally larger than the concentration of field 
halos of the same mass as a result of tidal mass-loss 
[55, 132, 133, 134, 135]. This is a small effect, of order 
10%, and it does not play a large role in the results we ob- 
tain below. This concentration-halo mass concentration 
has not just been established in numerical simulations: 
statistical weak lensing analysis of luminous red galaxies 
and clusters from the Sloan Digital Sky Survey (SDSS) 
has shown that the halo profiles of these objects are con- 
sistent with the NFW profile and mass-concentration re- 
lation given above [136]. 

The only remaining quantity in equation (30) that we 
haven't related to the mass of the subhalo yet, is the 
Sommerfeld enhancement factor S. The Sommerfeld en- 
hancement depends on the mass of the subhalo through 
the velocity dispersion. Therefore, we need to relate the 
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velocity dispersion of a DM subhalo to its mass. It is im- 
portant to note that all we really need is an approximate 
relationship between the velocity dispersion and the halo 
mass which fits well the range of Galactic size halos to a 
scale not much smaller than the mass scale of the smallest 
dwarf Spheroidals known today, since in most cases the 
Sommerfeld enhancement will saturate at velocity disper- 
sions of the order of the velocity dispersion of the dwarf 
Spheroidals. This means that the specific form that we 
adopt below for the velocity dispersion-mass scaling re- 
lation does not matter greatly. 

On the scale of clusters there is a well-defined rela- 
tionship between the mass of a DM halo and its one- 
dimensional velocity dispersion which one can derive 
from adiabatic scaling arguments [118] given by 

cr (X m^/^ . (34) 

This scaling relation can be easily derived in the case of 
an isothermal distribution function [137], since then 



p(r) 



(35) 



27rGr2 ' 

from which equation (34) immediately follows. This scal- 
ing relation has been found to describe numerical simu- 
lations of cluster-sized halos extremely well [77, 82, 118, 
138, 139, 140]. From recent numerical simulations it 
seems that we can use this relation even for subhalos of 
subhalos in a Galaxy sized dark matter halo. In Ref. [61] 
it is shown that the relation between the maximum cir- 
cular velocity of a subhalo and its mass is well fit by a 
power-law of the type 

m (X v^l^ . (36) 

Using the circular velocity at r_2 as a measure of the 
halo's velocity dispersion gives dy ex m^/^-^. Alterna- 
tively, we could use the circular velocity at the virial ra- 
dius as a measure of the velocity dispersion. In the case 
of an NFW profile we have that 

vcir^ir? m ' ^^^^ 

in which /(x) = ln(l + c)-c/(l + c) ^ 2.6(c/33)0-^ in the 
relevant concentration range. Using equation (33), this 
gives 

ay (X m°-^^ . (38) 

Therefore, we will use (as a rough approximation) the 
scaling relation given in equation (34). We calibrate this 
relationship at the mass scale of the dwarf spheroidals 
by using the velocity and mass of Canes Venatici I from 
Ref. [70]. 



D. The total and the "reduced" boost 

We define the total boost for a given subhalo mass as 
(see equation [29]) 



in which S{M) is the Sommerfeld enhancement for a halo 
of mass M and B{M) is the substructure boost for a 
halo of that mass, as defined in Eq. (30). This total 
boost for subhalos of mass M 10^ and M 10^ 
Mq calculated by following the procedure outlined in the 
previous subsection is shown in Figure 5 {top panels). 
It is immediately obvious that the total boost factors 
can get very large in certain regions of the parameter 
space spanned by m^/m^and a: boosts vary between 10^ 
to 10^. These masses correspond to the typical masses 
of the classical and newly-discovered dwarf Spheroidals 
(dSphs) and as such their velocity dispersions are of the 
order 10~^. The pattern of resonance regions and non- 
resonance regions that appears for the pure Sommerfeld 
enhancement (as shown in Figure 2) is also apparent in 
these plots of the total boosts, however, it is noteworthy 
that the resonances seems to be more highly peaked than 
they were in the pure Sommerfeld enhancement case, and 
that the boost the lower half of the parameter space, i.e., 
for smaller values of m^/m^, are higher than one would 
expect from the pure Sommerfeld enhancement combined 
with a boost factor from pure density enhancements. 

This is confirmed by looking at the "reduced boost" 
defined as 



B{M) 



S{M){l + Bo{M)) 



(40) 



B{M) = S{M) + B{M) , 



(39) 



in which Bo{M) is the boost we would find if we turned 
off the Sommerfeld enhancement in the calculation of the 
total boost, i.e., the boost we would get if we set S{m) 
equal to one in equation (30). This Bq corresponds to 
the boost factor that has been used before to estimate 
the total boost from density enhancements in dSphs and 
it is generally 0(10). This reduced boost therefore di- 
vides out the contribution of the Sommerfeld enhance- 
ment of the top-level subhalo and the boost from the 
over densities in lower- level halos, and it would be equal 
to one if there were no extra boost coming from the fact 
that the substructures are kinematically colder than the 
main subhalo. As we can see in Figure 5 {bottom panels), 
in various parts of the Yukawa coupling parameter space 
this reduced boost is equal to one such that there is no ex- 
tra boosts caused by the lower velocity dispersion in the 
substructures. This is mainly the case for large values of 
the parameter m^/m^, since the Sommerfeld enhance- 
ment, in the 1/v regime, levels off around S ~ rriff^/m^, 
and since the top-level subhalo has a velocity dispersion 
of ~ 10~^ this means that the Sommerfeld enhancement 
has mostly saturated at these large values of m^/m^, 
even in regions of parameter space in which a is small. 
Indeed, in the 1/v regime the Sommerfeld enhancement 
is ^ a/ j3 ~ 10^ for a = 10~^, such that the Sommerfeld 
enhancement has fully saturated over the whole a range 
for large values of m(f)/m^ {m(f)/m^ ~ 10~^). 

However, the Sommerfeld enhancement has not satu- 
rated in all regions of parameter space at this mass scale, 
and for smaller values of m^/m^ and especially in reso- 
nance regions, the reduced boost can reach high values. 
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FIG. 5: The total (defined in Eq. [39]) and reduced boost (defined in Eq. [40]) for subhalos of two different mass scales, 
dark matter profile of all the subhalos is assumed to be well-described by an NEW profile. 



The 



signaling a large extra boost caused by the kinematically 
coldness of the substructure. For small values of m^/m^ 
we see that we can get boosts that are an order of mag- 
nitude to two orders of magnitude larger than what we 
would expect from naively combining the Sommerfeld en- 
hancement from the velocity dispersion of the subhalo 
with the boost from density enhancements alone. This 
can be easily understood as being a consequence of the 
fact that the Sommerfeld enhancement has not saturated 
yet at the mass scale of the top-level subhalo, such that 
it can be boosted up to it saturation point by the contri- 
bution from the kinematically colder substructure. 

The largest values of the reduced boost are found in 
the resonance regions. In these regions the Sommerfeld 
enhancement both grows faster with decreasing relative 
velocity as well as saturates at a much larger value of the 
enhancement. Both of these effects conspire to give very 
large extra boosts. The resonance regions are clearly dis- 
tinguishable in this plot of the reduced boost as 45° lines. 



Because the resonance regions are so narrow the appar- 
ent extra boosts in different resonances seems to vary 
greatly, but this is mostly caused by the finite resolution 
of this figure. Extra boosts up to ~ 10^ are obtained in 
these resonances, and these extra boosts increase along 
the resonances when decreasing a, again caused by the 
higher saturation level for smaller m^/m^ values. 

The total boost factors obtained by assuming an NEW 
dark matter density profile do not depend much on the 
inner slope of the density profile. This is mainly because 
the calculation of the boost only considers relative lumi- 
nosities, for which the difference in normalization of the 
C oc plr^/D'^ relation for different inner slopes is irrele- 
vant, and because the mass of the subhalo is only weakly 
affected by the value of the inner slope — the function /(c) 
is only a weak function of the inner slope. This is con- 
firmed by calculating the total boost for a M = 10^ Mq 
subhalo for a subhalo with a very shallow inner cusp as 
well as for a subhalo with a very steep inner slope. These 
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FIG. 6: The total boost for a M = 10^ Mq subhalo for different shapes of the DM density profile: a generalized NFW profile 
described by an inner slope 7 {top panels); an Einasto profile described by a shape parameter a {bottom panels). 



are shown in Figure 6 {top panels), and, from comparing 
this with the total boost shown in Figure 5, it is clear that 
the inner slope does not affect the total boost caused by 
substructure, neither the overall magnitude of the effect 
nor the structure in the (m^/m^,^) plane. The main 
difference in DM annihilation flux between the different 
dark matter density profiles is therefore given by the nor- 
malization of the C oc plr^/D'^ relation, shown in Figure 
3. 



If the DM density profiles of halos is better described 
by an Einasto profile then the substructure boost is af- 
fected more strongly by the shape parameter a which 
basically describes how cored the DM density profile is. 
Changing a has a greater impact on changing the over- 
all DM density profile and therefore also plays a larger 
role in the relation between the characteristic density, 
the concentration, and the mass of the subhalo. The to- 
tal boosts for a M = 10^ Mq subhalo is shown in Figure 
6 {bottom panels) for two different values of the shape 
parameter a which span the range found in numerical 
simulations (for values of a greater than 0.2 the behavior 



basically saturates at that of a = 0.2). For a = 0.1 the 
total boost is about the same as for the GNFW profile 
shape, which is understandable given that the Einasto 
profile is more cuspy for smaller values of a. The same 
dependence on the parameter m^/m^ as is seen for the 
GNFW profile is displayed for the Einasto profile. In 
the case that a is larger the substructure boost reaches 
slightly larger values than for the GNFW density profile. 
Thus, as was the case for the GNFW profile, the impact 
of the different normalization of the 
ferent values of a will play a leading role in the difference 
between the different profile shapes. 

In Figure 7 the effect of varying two more of the many 
parameters introduced in the previous section is shown. 
The slope of the subhalo mass function n was set to - 
0.9 before, corresponding to the results from the latest 
high-resolution numerical simulations. In Figure 7 {top 
panels) we show the effect of both a shallower as well as 
a steeper subhalo mass function. The value of n = —0.7 
corresponds to the shallow slope found in Ref. [125]. We 
see that the effect of such a shallow slope is to seriously 
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FIG. 7: The total boost for a M = 10^ Mq for different slopes n for the subhalo mass function {top panels) and for different 
values of the dark matter thermal free streaming limit mo {bottom panels). 
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2.8 


5.2 X 10^ 


1,3 


Leo II 
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4 
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1.5 X 10^ 
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Sculptor 


79 
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6 
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2.8 
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Sextans 


86 
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2 


10 


1.0 X 10'^ 
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Ursa Major II 


32 
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2,4 


Ursa Minor 


66 
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TABLE I: Properties of the dwarf Spheroidals used in 


this si 


budy. 















References: (1) Penarrubia et al. [141], (2) Strigari et al. [67], (3) Mateo [69], (4) Simon and Geha [70]. 



reduce the substructure boost, which is simply a conse- ture is found in smaller and more concentrated clumps, 
quence of the fact that for a shallower slope less substruc- A steeper slope such as n = —1, as was found in the 
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FIG. 8: Predicted DM annihilation fluxes for various dwarf Spheroidal satellite galaxies. The background flux of extragalactic 
photons in the relevant energy range is given under the name of the dwarf Spheroidal. The properties of the dSphs used are 
given in Table I. 
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recent Via Lactea simulation [6] leads to a substructure 
boost which is about an order of magnitude larger than 
for the fiducial n = —0.9 value. 

The effect of the thermal free-streaming scale of DM 
on the substructure boost is shown in Figure 7 {bottom 
panels). This effect is relatively large since the free- 
streaming scale sets the scale of the smallest and most 
concentrated DM halos and raising or lowering this scale 
therefore changes an important part of the contribution 
to the substructure boost. For very large values of the 
free-streaming scale such as mo = 10~^ M©, the sub- 
structure boost becomes about an order of magnitude 
smaller. When the free- streaming scale becomes very 
small, much of the substructure is in very small, cold 
clumps, such that the total boost becomes about an or- 
der of magnitude larger. 



IV. PREDICTED ANNIHILATION SIGNALS 
FOR DWARF SPHEROIDALS 

Besides detection of DM annihilation from the Galac- 
tic Center or the diffuse component from the Galactic 
halo, there idea exists that the best candidates for DM 
annihilation detection might be the dwarf Spheroidal 
satellite galaxies (dSphs) of our Galaxy, since the dwarf 
galaxies are believed to be the most dark matter dom- 
inated structures in the Universe, and the DM annihi- 
lation signal coming from them would be much easier 
to distinguish from background sources than, e.g., would 
be the case for DM annihilation photons coming from 
the Galactic Center. That dSphs make good candidates 
for DM annihilation detection has been much discussed 
before, especially for the "classical", or pre-SDSS, dSphs 
[62, 110, 142, 143, 144, 145]. The annihilation signal com- 
ing from Draco has been the subject of much research in 
recent years because of its relative proximity to the Earth 
[65, 110, 144, 146, 147, 148, 149, 150]. The prospects 
for DM annihilation detection from dSphs has in general 
been found to be rather poor, with detection with Fermi 
being unlikely [150, 151]. Some of the dSphs have already 
been the subject of observational searches for DM anni- 
hilation and some limits on the DM annihilation signal 
have already been set [152, 153, 154, 155, 156, 157]. 

The dSphs discovered by SDSS have also been sug- 
gested to make good candidates for DM annihilation 
detection, because of their very large mass-to-light ra- 
tios (~ 10^) [67, 68]. In the context of the Sommerfeld 
enhancement, these recently discovered dSphs are espe- 
cially interesting, since they also have the smallest veloc- 
ity dispersion of all of the known satellites of our Galaxy. 
Predicting the DM annihilation signal from these low- 
surface brightness galaxies is complicated by the dearth 
of kinematical tracers of the DM distribution and their 
contamination with background stars [158]. Neverthe- 
less, detailed studies have constrained the DM density 
and velocity dispersion profiles of some of these dSphs 
[65, 73, 141]. We will use some of the results on the 



DM density and velocity profiles found in these studies 
to show the extra effect of the Sommerfeld enhancement 
on the predicted DM annihilation signal from some of 
these dSphs. 

We consider the dSphs given in Table I here. Eight of 
the dSphs are classical dwarf galaxies and we use the 
DM density profiles fitted to them in Ref. [141]. In 
this the DM density profiles of the dSphs were not com- 
pletely determined by the available data — a degeneracy 
between the parameters of the NEW profile used to de- 
scribe the DM density and the King profile used to fit the 
luminous component of the galaxies was found — and the 
concentration-mass relation from CDM [71, 72] was used 
to break the degeneracy. Since these authors expressed 
their results in the form of the rmax and ^max param- 
eters of the DM density profiles, we use the procedure 
outlined in Section IIIB to convert these to ps and Ts 
values. As these 'i^max and r^iax parameters were fit to 
an NEW DM profile, we will only consider an NEW pro- 
file for the dSphs here, but the general conclusions about 
the importance of the exact profile shape also hold in 
these particular cases. The velocity dispersions of these 
dSphs are taken from Ref. [69] and we again model these 
subhalos as having a constant velocity dispersion. 

We also consider some of the newly discovered dSphs. 
In particular we discuss the prospects for DM annihila- 
tion detection from Coma Berenices and Ursa Major II, 
and we use the best fit values to the available data from 
Ref. [67]. Since these authors directly quote ps and rg 
values we do not include r^iax and ^max values for these 
dSphs. We use the velocity dispersions determined for 
these galaxies in Ref. [70]. 

In order to calculate the substructure boost to the DM 
annihilation cross section we need the mass of the dSphs. 
Many definitions of the mass of the galaxy could be used 
here, among the more appropriate ones is the mass con- 
tained inside the tidal radius, i.e., the radius inside of 
which material in the dSph is safe from being tidally 
stripped by the Galaxy. However, since we are consid- 
ering the DM annihilation flux coming from within the 
radius r_2 we will use the mass contained within this ra- 
dius to calculate the substructure boost. This is more ap- 
propriate since (1) it can be calculated straightforwardly 
from the DM density profile of the dSph and (2) the mass 
contained within r_2 is really the more relevant mass to 
calculate the substructure boost, since we are only inter- 
ested in contributions to the DM annihilation flux coming 
from within the radius r_2. 

Because of the similarities between some of the DM 
density profiles of the dSphs in Table I we do not show 
the predicted annihilation fluxes for all of the dSphs in 
Figure 8, instead focusing on a representative sample of 
them. Inspection of Table I shows that annihilation flux 
of the dSphs not shown here will be similar to that of the 
ones shown in Figure 8. 

At this point we can also make a connection to another 
quantity that is often used to describe the astrophysi- 
cal contribution to the DM annihilation signal, J. This 
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quantity is generally defined as 



J = 



J An J\os 



(41) 



and is therefore equal to the hne-of-sight integral over 
the density squared averaged over solid angle. We can 
define a slight variant of this which takes into account 
the Sommerfeld enhancement 



J 



1 

An 



An J\os 



(42) 



in which S is the Sommerfeld enhancement averaged over 
the distribution of relative velocities. If we work this 
out for the dwarf Spheroidals, using the expression from 
equation (14) for the line-of-sight integral averaged over 
the solid angle spanned by the scale radius we find that 
this J is approximately equal to the total boost S, with 
a pre-factor 0{1). Therefore in the case of the dSphs 
studied here the J can be approximately read off from 
the top panels of Figure 5. 

It is clear from Figure 8 that the DM annihilation fiux 
predicted to be coming from the dSphs is very large in 
a large class of models. We have to stress here that our 
predicted fiuxes are all for the particle physics specific pa- 
rameters given in Section III A, i.e., (av) = 3 x 10~^^ cm^ 
s~^, = 700 GeV and a particular form for the DM 
annihilation spectrum. To predict the gamma-ray fiux 
for a particular model with different values of these par- 
ticle physics parameters one is required to calculate the 
particle physics parameter specific factor to the DM an- 
nihilation fiux and appropriately rescale the values given 
here. In order to compare these signals to the expected 
background fiux we will only consider contributions to 
the background fiux of extragalactic gamma rays. Al- 
though there are also contributions to the background 
fiux from pion-decay, inverse Compton scattering, and 
Bremsstrahlung, at the high Galactic latitudes at which 
the dwarf Spheroidals are found these are all of the same 
order or smaller than the extragalactic background [159]. 
We have integrated the extragalactic gamma ray fiux as 
measured by EGRET over the energy range between 4 
GeV to 250 GeV and multiplied it with the solid an- 
gle that the dSphs span on the sky. The extragalactic 
gamma ray fiux is given by a power law [160] 



dN^ _ f E_ 
dEdAdtdVt ~ \E~Q 



(43) 



where /c = 7.3 x 10~^ cm~^ s~^ sr~^ GeV~^, a = 
2.10, and Eq = 0.45 GeV. We have included the back- 
ground fiux of extragalactic photons in Figure 8. This 
background fiux is in all cases of the order of 10~^ to 
10~^^ photons cm~^ s~^, such that the DM annihilation 
boosted by substructure and Sommerfeld enhancement 
outshine this background in most of the parameter space, 
and by a few orders of magnitude in most models. Given 
that the fiux above 100 MeV over a 55 days time-span 
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FIG. 9: The distances to and masses of subhalos within the 
tidal radius of the main halo in Via Lactea II. The distances 
are for three random observers at a distance of 8.5 kpc from 
the galactic center of the main halo. 



needed to separate a point source of DM annihilation 
from the background by Fermi to give a 5cr detection 
is about 10~^ photons cm~^ s~^ and that the necessary 
fiuxes scale roughly as texp [114], the DM annihilation 
of these dSphs should be detectable by Fermi in many 
of the models parametrized by (m^/m^,^). The Ursa 
Major II dSph is the only dSph for which the DM anni- 
hilation dominates the background in all the regions of 
parameter space. This shows that the newly discovered, 
most-dark-matter-dominated satellite galaxies are possi- 
bly the best targets for the indirect detection of DM. We 
conclude that the prospects for detecting DM annihila- 
tion from the dSphs with Fermi are very good in models 
that include a significant Sommerfeld enhancement. 



V. THE ANNIHILATION SIGNAL FROM 
SUBHALOS IN VIA LACTEA II 

Besides the satellite galaxies that have been observed 
in the Galactic halo CDM predicts that the halo is filled 
with hundreds of individual subhalos — a fact which gives 
rise to the so-called "missing satellites problem" [57, 58]. 
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FIG. 11: Total gamma-ray flux from subhalos in Via Lactea 
II for 1,000 random observers at 8.5 kpc from the center of the 
main Via Lactea halo for m^/m^ = 10~^ and a = 10""^. For 
comparison we show the total gamma-ray flux from subha- 
los in the absence of Sommerfeld enhancement as the dashed 
histogram and the top horizontal axis. 



FIG. 10: Top: Histogram of the predicted gamma-ray flux 
coming from individual subhalos contained within the tidal 
radius of the main halo in Via Lactea II for three random 
observer placed at a distance of 8.5 kpc from the center of 
the main Via Lactea halo for mcj^/m^ = 10~^ and a = 10~^. 
Bottom: Histogram of the total boosts of the subhalos for the 
same values of m(f)/m^ and a. 



Each of these subhalos has a mass of the order of the 
masses of the observed dSphs or smaller, and thus will be 
kinematically cold with respect to the main halo. Thus, 
the total DM annihilation signal coming from all of the 
substructure will be boosted by the Sommerfeld enhance- 
ment described above, such that an appreciable DM an- 
nihilation signal could be expected from all of the sub- 
structure. In this section we get a feeling for the expected 
signal from DM substructures in a typical Galaxy-sized 
halo. If the DM annihilation signal from subhalos with- 
out star formation could be detected this could give a 
clear signal that the DM exists and would open up a 
whole new arena for testing CDM. 

We use the distribution of DM subhalos from the Via 
Lactea II numerical A^-body simulation^ [38] . Via Lactea 
II is a state-of-the-art TV-body simulation of a Galaxy- 
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sized halo since redshift 104.3. It sampled the region 
of the main halo with 1.1 x 10^ particles, each with a 
mass of 4,100 Mq. Over 40,000 subhalos can be resolved 
within 402 kpc of the center of the main halo and 20,047 
subhalos are detected which had a peak circular velocity 
larger than 4 km s~^ at some time. The properties of 
the subhalos in this latter sample can be accurately es- 
tablished, i.e., they are unaffected by resolution effects. 
We select from this sample the 9,830 subhalos which lie 
within the tidal radius of the main halo, which is ap- 
proximately 462 kpc. We show the distribution of these 
subhalos in mass and distance from the Sun in Figure 9. 
We have placed the Sun at random locations at a distance 
of 8.5 kpc of the center of the main halo, corresponding 
to an earth-based-observer-like vantage point. The dis- 
tribution of the distances in Figure 9 is given for three 
of these random observers and the distribution is very 
similar for all three observers. We see that most of the 
subhalos have masses in the range 10"^ to 10^ Mo, with 
subhalos with masses of the order of the dSphs being rel- 
atively rare. The lower limit is of course set by the finite 
numerical resolution of the simulation and doesn't corre- 
spond to a physical cut-off — indeed, the physical cut-off 
lies at the thermal free-streaming scale. The distribution 
of the distance to the subhalos shows that the distribu- 
tion of the subhalos is strongly radially anti-biased, with 
a sharp cut-off caused by our cut at the tidal radius of 
the main halo. This radial anti-bias is a consequence of 
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FIG. 12: Same as Figure 10 but for mf^jm^ = 10 ^ and FIG. 14: Same as Figure 10 but for a resonance: mf^jr. 
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the fact that there is much more space far away from the 
center of the galaxy, and because of the strong tidal forces 
near the center of the galaxy, which strongly affects halos 
which pass near the center of the galaxy. 

For each halo we can estimate the velocity dispersion 
using the mass- velocity dispersion relation given in equa- 
tion (34) and we can calculate the characteristic density 
ps and the characteristic length-scale from the values 
of ^max ciiid Tmax given for each subhalo, by the proce- 
dure outlined at the end of Sec. IIIB. For a given set 
of particle physics parameters m^/m^ and a this allows 
us to calculate the total boost to the DM annihilation 
cross section for each of the subhalos, and by using the 
assumptions about the particle physics given in Sec. Ill A 
we can also calculate the photon flux from DM annihila- 
tion coming from each subhalo, and the total gamma-ray 
flux coming from all of the subhalos combined. 

We only cover some typical cases for the Sommer- 
feld enhancement here which suffices to show the gen- 
eral trend. We consider the particle physics parameters 
rriff) / and a which give rise to the smallest Sommerfeld 
enhancements in Figures 10 and 11, i.e., mcjy/m^ = 10~^ 
and a = 0.001. The top panel of Figure 10 shows the 
distribution of gamma-ray fluxes from DM annihilation 
for the subhalos, while the bottom panel shows the dis- 
tribution of the total boost factors. As can be seen, the 
total boost factors are quite small: the typical boost fac- 
tor is only about 30 with the highest boost factors being 
about 100. The total gamma-ray flux coming from sub- 
halos is then also quite small in this case, as is shown in 
Figure 11, in which the total gamma-ray flux is shown for 
1,000 random observers at a distance of 8.5 kpc from the 
galactic center. At most observation locations this total 
flux is about lO"*" cm^ s~^, with the high-end tail not 
extending further than about half an order of magnitude 
from this. For comparison the dashed histogram in Fig- 
ure 11 shows the same distribution of total gamma-ray 
flux for different random observers in the absence of any 
Sommerfeld enhancement. This confirms that the total 
gamma-ray flux with Sommerfeld enhancement is about 
an order of magnitude larger than without Sommerfeld 
enhancement. 

Figures 12 and 13 show the same for a different ex- 
tremum of the Sommerfeld enhancement: the case in 
which m^/m^ = 10~^ and a = 0.1. As we saw before, 
in this case the Sommerfeld enhancement can reach very 
high values, and thus we see that the typical DM an- 
nihilation signal coming from DM subhalos is large: the 
typical flux for one subhalo is equal to the total flux in the 
scenario described in the previous paragraph. The total 
boost factors for most of the halos are very large, of the 
order of 10^, which gives rise to this very large photon- 
flux. The total predicted flux in this case is many orders 
of magnitude larger than the total flux in the other ex- 
tremum of the Sommerfeld enhancement and about seven 
orders of magnitude larger than the total gamma-ray flux 
in the absence of any Sommerfeld enhancement. Thus, 
the exact value of the parameters m^/m^ and a, upon 



which the Sommerfeld enhancement depends, matters a 
great deal for the predicted gamma-ray flux. 

Finally, we look at the expected signal coming from 
subhalos when the Sommerfeld enhancement is at a reso- 
nance, the same distribution as before are shown in Fig- 
ures 14 and 15. Again we find large boost factors and 
large individual and total gamma-ray fluxes, although 
they are smaller than the values found for the / = 
10 ~^ and a = 0.1 parameter set. However, since these 
values are near a resonance in a region of parameter space 
around which the predicted signal looks more like that 
for the parameter set m^/m^ = 10~^ and a = 0.001, this 
again shows that whether dark halos, i.e., subhalos with- 
out star formation, can be detected through DM annihi- 
lation products depends very sensitively on the particle 
physics parameters upon which the Sommerfeld enhance- 
ment depends. 

In conclusion we can say that in a scenario in which 
there is Sommerfeld enhancement there is much more 
hope to detect dark subhalos than in more conventional 
scenarios, but that the details of the prediction depend 
strongly on the underlying particle physics. 



VI. CONCLUSION 

The recently proposed Sommerfeld enhancement to the 
DM annihilation cross section forces us to reconsider 
statements made in the past about the expected signal of 
DM annihilation from cold substructures in the halo of 
the Galaxy. Since the Sommerfeld enhancement depends 
critically on the velocity dispersion of the substructure in 
question, as well as on the details of the particle physics 
model used to describe the attractive (or, possibly, re- 
pulsive) force that gives rise to the Sommerfeld enhance- 
ment, the contributions from the specific particle physics 
model and the details of the DM distribution in the sub- 
structure can no longer be as cleanly separated as they 
have been before. This means that in a scenario in which 
there is a significant Sommerfeld enhancement, previ- 
ous claims of "particle physics independent" conclusions 
about the DM annihilation signal from substructure (e.g., 
[39, 65]) need to be reconsidered in this new framework. 
This complicates the analysis but it opens up new possi- 
bilities for detecting DM annihilation and learning about 
the properties of the dark sector. As we showed above, 
the combination of the Sommerfeld enhancement and the 
substructure that is believed to exist down to small scales 
leads to very large boost factors because of the combined 
effect of the higher densities of small subhalos as well 
as the fact that they are kinematically colder than the 
larger halos. This leads to predictions for the DM anni- 
hilation gamma-ray flux that are much larger than was 
previously believed. As we showed in Sec. IV, in mod- 
els in which there is a significant Sommerfeld enhance- 
ment at the level of the Galactic halo and its subhalos, 
the DM annihilation signal from the dwarf Spheroidals 
outshines the extragalactic background flux and the pre- 



19 



dieted gamma-ray fluxes are so large that they should be 
detectable by Fermi. Given a detailed knowledge of the 
density structure and velocity dispersion profile of the 
dwarf Spheroidal galaxies, much can then be learned if 
we are able to detect the DM annihilation signal com- 
ing from these dSphs, by considering the relative signal 
and using the results of Sec. IV, or by using the detailed 
velocity dispersion profile of the dSphs [161]. 

We have focused on the consequences of the Sommer- 
feld enhancement for DM subhalos in a ACDM cosmol- 
ogy, using the latest results from high-resolution numeri- 
cal A^-body simulations to inform our assumptions about 
the abundance and properties of the substructure. Un- 
der these assumptions we were able to show that there 
are very large boosts to the DM annihilation cross sec- 
tion from the combined effect of substructure boosts and 
Sommerfeld enhancement in large parts of the parame- 
ter space consisting of the parameters of the attractive 
Yukawa type force between the DM particles which is 
responsible for the Sommerfeld enhancement. These re- 
gions of parameter space are in no way unrealistic. In- 
deed, in one of the most straightforward interpretations 
of the ATIC/PAMELA results, the sharp cut-off around 
1 TeV seen in the ATIC spectrum is interpreted as being 
a consequence of the fact that the dark matter mass is 
about 700 GeV. The absence of an excess in anti-protons 
similar to the excess in positrons seen by PAMELA can 
be explained by positing that the dark matter particles 
annihilate to an intermediate boson (/), which is lighter 
than a proton (or even, lighter than the effective quark 
mass Aqcd) in order to kinematically forbid the an- 
nihilation into protons [30]. This naturally leads to a 
m(f)/m^< 10~^. The coupling constant of the new force 
in the dark sector is also generally a > 10~^ [8, 26]. Thus, 
from Figure 5, we expect there to be a rather large boost 
to the DM annihilation cross section in these models. 
This leads to particularly large DM annihilation signals 
predicted for the dSphs (see Figure 8). 

We have not discussed the substructure boost to the 



DM annihilation signal from unresolved substructure, 
such as very small subhalos and DM tidal debris from 
disrupted subhalos, in the Galactic halo. This is, how- 
ever, interesting in the light of the Sommerfeld enhance- 
ment, since it might lead to a different prediction about 
the relative magnitudes of the annihilation signal com- 
ing from the resolved Galactic halo and unresolved halo, 
as well from the resolved substructures in the Galactic 
halos. It is already clear that in the generic 1/v Som- 
merfeld enhancement the main halo, with its typical ve- 
locity dispersion of ~ 150 km s~^, will receive a much 
smaller enhancement than the dSphs, with their typical 
velocity dispersions of 10 km s~^. This, however, ignores 
the effects of saturation and resonances, which might be 
important in the specific particle physics model adopted 
for the Sommerfeld enhancement. Additionally, while 
the velocity dispersion of the Galactic halo might be of 
the order of 100 km s~^, the actual velocity dispersion 
relevant to the Sommerfeld enhancement might be much 
smaller for halo DM particles. If a large fraction of the 
DM in the Galactic halo is part of a cold stream — and 
various numerical simulations predict that large parts of 
the halo, if not all of the halo, are made up of such cold 
streams with typical velocity dispersions of the order of 
1 km s~^ [59, 60] — then, in addition to any boost caused 
by density enhancements in the stream [162], the boost 
from the Sommerfeld enhancement could also play a large 
role, as the relevant relative velocity of DM particles is 
set by the velocity dispersion of particles in the stream, 
and not by the velocity dispersion of the streams in the 
halo, which is what the 100 km s~^ Galactic halo veloc- 
ity dispersion corresponds to in these scenarios. Thus, 
the expected DM annihilation signal from the main halo 
could be much larger than that naively expected for a 
halo with a velocity dispersion of 100 km s~^. 
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